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PURPOSE:  This  Coastal  and  Hydraulics  Laboratory  technical  note  describes  the  process  of 
scaling  sediment  for  the  1 :20  General  Physical  Model  of  the  Fish  Barrier  Structure  downstream 
of  Mud  Mountain  Dam,  which  is  located  on  the  White  River  near  Seattle,  WA. 

INTRODUCTION:  The  White  River  is  a  mountainous  river  with  head  waters  on  Mount  Rainier. 
The  Fish  Barrier  Structure  on  the  White  River  serves  as  a  collection  point  of  migratory  fish. 
Operators  at  the  structure  collect  the  fish  and  transport  them  upstream  of  Mud  Mountain  Dam. 
Currently,  the  fish  barrier  uses  a  weir  to  create  a  head  differential,  which  stops  the  passage  of  the 
migrating  fish.  However,  the  head  difference  causes  additional  deposition  upstream  of  the  fish 
barrier.  Furthermore,  with  large  annual  sediment  loads,  the  White  River  is  naturally  aggrading, 
and  the  fixed  weir  does  not  allow  for  structure  operations  to  mobilize  the  sediment.  Thus,  a  new 
structure  is  desired  to  both  mitigate  sediment  accumulation  and  increase  fish  collection  capacity. 

Correctly  selecting  the  appropriate  physical  model  sediment  is  critical  for  properly  simulating  the 
prototype  behavior  of  the  proposed  structure  (Franco  1978).  Numerous  studies  (Carr  et  al.  2015; 
Einstein  and  Chien  1956;  Henriquex  et  al.  2008;  Kocyigit  et  al.  2005;  Lu  et  al.  2013;  Yalin  1989) 
have  expended  extensive  effort  in  an  attempt  to  replicate  prototype  sediment  behavior  in 
movable  bed  models.  Properly  capturing  the  movable  bed  freedom  is  imperative  when  channel 
modifications  or  complex  flow  patterns  alter  bed  alignment.  This  freedom  is  directly  dependent 
on  appropriately  scaled  model  sediment. 

For  this  model,  there  are  two  basic  sediment  behaviors  that  are  of  principal  interest: 

•  the  stability  of  the  coarse  sediment  bed  (i.e.,  probability  of  erosion/scour) 

•  the  deposition  of  fine  sediment  from  suspension  (wash  load). 

To  represent  these  two  behaviors  in  a  single  physical  model,  it  was  determined  that  the  best 
approach  was  to  select  two  separate  sediment  classes  to  be  modeled:  one  bed  material  class  that 
is  scaled  for  initiation  of  motion  and  one  suspended  sediment  class  that  is  scaled  for  deposition 
from  suspension. 

EFFECTIVE  GRAIN  SIZE:  The  inherent  complexities  associated  with  sediment  gradations 
cannot  (in  general)  be  properly  scaled  in  a  physical  model.  Therefore,  a  simplification  is  introduced 
whereby  an  effective  grain  size  is  estimated  for  the  prototype  mixture.  The  effective  grain  size  is 
the  size  required  for  a  single  grain-sized  model  sediment  to  behave  the  same  as  a  prototype  sedi¬ 
ment  gradation,  with  respect  to  some  desired  characteristic  (Christensen  1969).  For  the  purposes  of 
this  model,  the  desired  characteristic  is  the  probability  of  motion.  Once  an  effective  grain  size  is 
defined  based  on  the  prototype  gradation,  this  effective  grain  size  is  scaled  for  model  application. 
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The  sediment  scaling  analysis  for  this  effort  was  guided  by  the  need  to  determine  which 
prototype  sediment  classes  would  be  modeled  as  bedload  and  which  would  be  modeled  as 
suspended  load.  The  prototype  sediment  is  generally  subdivided  into  two  groups:  wash  load 
(WL)  and  bed  material  load  (BML).  These  two  components  comprise  the  total  sediment  load. 
The  WL  is  the  sediment  in  transport  derived  from  sources  other  than  the  bed  (Lane  1947).  The 
WL  is  finer  than  the  BML  and  is  not  found  in  appreciable  quantities  in  the  bed  (Biedenham  et  al. 
2006).  Conversely,  the  BML  is  sediment  whose  particles  are  found  in  appreciable  quantities  in 
the  bed  (Vanoni  2006;  Biedenham  et  al.  2006).  However,  the  delineation  of  what  percentage 
constitutes  an  appreciable  quantity  is  abstract  and  challenging  (Biedenham  et  al.  2006).  Because 
the  definitions  of  WL  and  BML  are  qualitative,  additional  interpretation  to  discretely  define 
these  components  is  required. 

For  this  effort,  the  delineation  between  WL  and  BML  is  defined  by  the  criteria  of  Einstein  (1950) 
as  the  D10  (grain  size  in  which  10%  of  the  particles  are  smaller)  of  the  mixture,  in  which  WL  and 
BML  are  finer  and  coarser  than  10%  of  the  sediments  by  gradation,  respectively.  This  analysis 
should  be  based  on  a  representative  prototype  gradation  taken  from  non-slack  water  areas  (Einstein 
1950).  For  this  study,  the  sediment  sampling  was  limited  to  locations  downstream  of  the  existing 
fish  barrier  structure.  Figure  1  presents  the  sieve  analysis  for  two  downstream  samples  (500  feet 
downstream  of  barrier)  of  the  prototype  sediment.  Taking  the  (numerical)  derivative  of  the  percent 
passing  curve  of  the  representative  prototype  gradation  results  in  the  fraction  of  each  class  in  the 
sample  by  sieve  size.  The  fractional  analysis  reveals  a  clear  bimodal  distribution  of  the  prototype 
sediment  (Figure  2).  This  bimodal  distribution  suggests  that  simulating  the  sediment  behavior  of 
the  prototype  with  two  separate  sediment  classes  is  appropriate. 


Figure  1 .  Representative  sediment  sample  (eight  samples)  collected  from 
downstream  of  the  existing  Fish  Blockage  Structure,  White  River 
WA. 
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Figure  2.  Percent  retained  from  representative  sediment  sample. 


The  bimodal  break  was  2  millimeters  (mm),  which  corresponds  to  the  DIO  of  one  of  the  samples 
(Figure  1).  Here,  the  coarser  material  (>  DIO)  is  designated  as  the  BML  while  the  finer  material 
(<  DIO)  is  designated  as  the  WL.  By  separating  the  two  quantities  into  their  own  respective 
gradations  and  using  Equation  1  (Christensen  1969),  the  effective  grain  size  (de)  is  calculated  for 
each  representative  class.  The  effective  grain  sizes  for  the  White  River  downstream  of  the  Mud 
Mountain  Fish  Barrier  Structure  are  found  from  eight  samples  to  be  29.5  mm  (coarse  gravel)  and 
0.27  mm  (medium  sand)  for  BML  and  WL,  respectively  (see  nomenclature  at  end  of  text  for 
variable  definitions). 


SIMILITUDE:  Free  surface  flow  similitude  can  be  defined  as  follows: 


We  say  that  two  flows  are  similar  if  they  are  geometrically  similar  and  if  all  the  relevant 
dimensionless  parameters  are  the  same  for  both  flows.  (Stockstill  2015) 

Similitude  principles  form  the  design  fundamentals  such  that  the  results  of  the  model  can  be 
converted  to  prototype.  For  any  model,  free  surface  similitude  requires  the  existence  of 
geometric,  kinematic,  and  dynamic  similitude.  Geometric  similitude  exists  between  the  model 
and  prototype  if  the  ratios  of  all  corresponding  dimensions  in  the  model  and  prototype  are  equal. 
Kinematic  similitude  exists  between  model  and  prototype  if  the  flow  paths  of  the  moving 
particles  are  geometrically  similar  and  if  the  ratios  of  the  particle  velocities  are  equal.  Dynamic 
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similitude  exists  within  geometrically  and  kinematically  similar  systems  when  the  ratio  of  all 
homologous  forces  in  the  model  and  prototype  are  equal. 

Unfortunately,  it  is  not  possible  to  achieve  true  dynamic  similitude  for  all  forces  in  a  scaled 
physical  model  (ASCE  2000).  However,  most  physical  systems  are  dominated  by  a  small 
number  of  the  forces.  If  two  dominant  forces  can  be  identified,  the  non-dimensional  ratio  of 
these  forces  can  be  scaled,  ensuring  dynamic  similitude  of  these  forces. 1 

The  non-dimensional  term  that  defines  dynamic  similitude  for  free-surface  flow  is  the  Froude 
number,  where  the  ratio  of  inertia  to  gravity  is  equal  between  the  model  and  prototype.  Particle 
similitude  can  be  evaluated  in  terms  of  two  governing  force  ratios: 

•  incipient  motion  (Shields  similitude) 

•  suspension  and  settling  (Rouse  similitude)  (ASCE  2000). 

Julien  (2002)  outlined  four  important  parameters  for  moveable  bed  similitude: 

•  Froude  similitude 

•  resistance/roughness  similitude;  bed  roughness  similitude  is  critical  for  similar 
hydrodynamics  from  prototype  to  model  such  as  water  slope  and  velocity 

•  similitude  of  dimensionless  grain  diameter,  defining  sediment  grain  regime 

•  similitude  of  the  bed  entrainment  of  individual  grains  into  suspension. 

By  correctly  accounting  for  these  four  parameters  defined  by  Julien  (2002),  it  is  possible  to 
acquire  physical  model  results  applicable  to  the  prototype.  However,  perfect  similitude  for 
sediment  and  flow  is  unattainable  in  a  physical  model,  so  the  dominant  forces  must  be  identified 
and  properly  scaled  (i.e.,  enforced).  For  this  study,  Shields  similitude  will  be  enforced  for  the 
BML  (to  ensure  that  the  initiation  of  motion  is  scaled  properly),  and  Rouse  similitude  will  be 
enforced  for  the  WL  (to  ensure  that  suspension  and  deposition  are  scaled  properly). 

The  selection  of  model  sediment  classes  that  satisfy  these  criteria  requires  the  use  of  four 
dimensionless  numbers: 

•  Rouse  number  (or  particle  Peclet  number)  (Equation  2) 

•  Shields  parameter  (Equation  3) 

•  dimensionless  grain  diameter  (Equation  4) 

•  particle  Reynolds  number  (Equation  5). 

The  equations  describing  these  terms  are  given  below.  The  following  sections  will  discuss  how 
they  are  used  to  select  the  model  sediments. 


P 


(2) 


1  Note  that,  for  any  physical  model,  the  neglected  forces  are  still  important  as  they  determine  the  physical  limitations  on  scale. 
For  example,  Froude  models  must  be  large  enough  to  remain  turbulent  (Reynolds  constraint)  and  to  ensure  that  surface 
tension  does  not  become  significant  in  the  model  (Weber  constraint). 
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BED  MATERIAL  LOAD  SCALING:  By  applying  the  effective  grain  size,  the  BML  for  the 
prototype  was  scaled  for  initiation  of  motion.  The  defining  non-dimensional  parameter  for 
initiation  of  motion  is  the  Shields  parameter  (Shields  1936)  (Figure  3).  It  is  common  practice  in 
physical  modeling  to  assume  that  the  Shields  parameter  is  a  constant  between  model  and 
prototype:  that  is,  the  critical  shear  stress  in  the  model  and  the  prototype  scales  only  with  the 
grain  diameter  and  with  the  specific  gravity  of  the  sediment.  For  this  study,  the  prototype 
sediment  consists  of  gravel  and  cobble.  Geometrically  scaling  the  prototype  effective  grain  size 
results  in  a  sand  size  gradation.  A  sand  and  gravel-based  sediment  does  not  in  general  have  the 
same  Shields  parameter.  Therefore,  it  is  necessary  to  ascertain  the  appropriate  Shields  parameter 
for  the  prototype  sediments. 


The  most  direct  way  to  obtain  the  Shields  parameter  is  to  determine  directly  the  critical  shear 
stress  of  both  the  prototype  sediments  (BML  and  WL)  by  performing  flume  tests  and  then 
calculate  the  Shields  parameter  from  Equation  3.  However,  for  this  study,  flume  testing  of  the 
prototype  sediment  mixture  was  not  a  viable  option  since  the  D90  of  the  prototype  is  nearly  180 
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mm  and  would  require  an  estimated  unit  discharge  35  cubic  feet/second/foot  to  mobilize.  This 
discharge  is  beyond  the  discharge  capacity  of  the  available  test  facilities. 

Since  direct  measuring  of  the  critical  shear  stress  of  the  prototype  sediment  proved  impractical,  a 
novel,  alternative  method  for  estimating  the  Shields  parameter  was  developed  for  this  study.  It  can 
be  shown  from  first  principles  that  the  Shields  parameter  is  a  function  of  the  angle  of  repose  of  the 
sediment  (Christensen  1993).  The  angle  of  repose  is  defined  as  the  maximum  or  critical  angle 
(Bagnold  1941)  when  a  granular  mass  destabilizes  and  slides  (Carrigy  1970;  Deganutti  et  al.  201 1). 
Several  methods  for  determining  the  Shields  parameter  indicate  the  dependency  of  the  Shields 
parameter  on  the  angle  of  repose  (e.g.,  Julien  [1995]).  Therefore,  for  a  given  sediment  mixture,  the 
Shields  parameter  can  be  estimated  indirectly  by  measuring  the  angle  of  repose  for  the  mixture. 

To  measure  the  angle  of  repose  of  the  prototype  sediment,  a  tilt- tank  apparatus  was  designed  and 
fabricated.  The  tilt-tank  testing  apparatus  is  an  acrylic  box  with  a  removable  and  sealable  top. 
The  sediment  sample  is  loosely  placed  in  the  box  and  submerged  in  water.  Once  the  sample  is 
saturated,  the  box  is  tilted  at  a  rate  of  0.2  degrees  (deg)/second.  Particle  sifting  can  occur,  but  this 
is  not  considered  failure.  A  failure  is  said  to  occur  when  multiple  individual  particles  are 
dislodged  at  the  same  time.  Table  1  shows  the  synopsis  of  the  results  for  a  series  of  1 1  tests  of 
the  prototype  sediment  conducted  for  this  study. 


Table  1.  Tilt-tank  test  results  for  the  prototype  bed  material  load. 

Prototype  Bed  Material  Load  Analysis 

Bed  Layer 
Thickness,  mm 

Average  Angle  of 
Repose,  deg 

Standard 
Deviation,  deg 

Gradation,  mm 

D10 

|  D50 

|  D90 

165.1 

46.4  ] 

2.49 

0.7 

27.9 

125.6 

241.3 

45.3  ] 

3.55 

0.6 

27.9 

113.4 

The  values  for  the  angle  of  repose  are  larger  than  those  reported  by  Deganutti  (2011).  For  a 
mixture  of  22.4  <  dso  <  31.5  mm,  the  average  angle  of  repose  was  40.8  deg  (Deganutti  2011). 
However,  the  gradation  was  more  uniform  and  was  described  as  rounded.  Carrigy  (1970)  showed 
that  the  angle  of  repose  for  a  uniformed  submerged  natural  material  was  less  than  the  same 
material  in  air  (also  noted  by  Land  [1964]).  These  gradations  were  not  as  widely  distributed  as 
the  present  prototype  material. 

Once  the  angle  of  repose  was  defined  for  the  prototype  material,  the  model  material  was  scaled. 
First,  the  Shields  number  for  the  prototype  material  was  calculated  using  Equation  3.  Likewise,  it 
was  calculated  for  the  model.  However,  an  initial  guess  at  the  effective  grain  size  of  the  model  was 
required.  This  was  done  by  geometrically  scaling  the  prototype  effective  grain  size  by  X  (1:20). 
Once  both  the  model  and  prototype  Shields  number  were  estimated,  Equation  6  was  applied: 


d 


m 


(6) 


This  provided  a  new  estimate  for  dm  (model  sediment  diameter)  and  was  substituted  into 
Equation  4.  The  dimensionless  grain  size  was  calculated  and  used  in  Equation  7  Julien  (2010): 
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T*cm  =  0.06ton(J)m  When  dt  >  50  (7) 

Then  the  new  estimate  for  the  models  Shields  number  was  re-applied  in  Equation  6  to  repeat  the 
iteration.  The  iteration  was  looped  until  dm  stopped  changing. 

In  order  to  characterize  the  sensitivity  of  the  above  calculation  to  uncertainties  in  the  input 
parameters,  a  perturbation  analysis  was  performed.  The  variation  in  particle  sizing  that  resulted 
from  this  perturbation  analysis  is  shown  in  Table  2.  Perturbations  of  both  the  effective  grain  size 
and  angle  of  repose  of  the  prototype  sediments  were  performed.  The  effective  grain  size  was 
varied  by  +/-  20%  and  +/-  40%.  Justification  for  the  wide  range  in  variation  in  effective  grain 
size  arises  from  the  observed  tendency  for  skewing  of  large  gradation  sizes  toward  the  smaller 
end.  This  skew  is  due  primarily  to  the  large  granular  material  and  the  difficulty  in  obtaining  a 
representative  sample.  Variation  in  the  angle  of  repose  was  based  on  the  mean  (45.3  deg)  and  +/- 
one  standard  deviation  found  from  the  tilt-tank  testing. 


Table  2.  Sensitivity  of  scaled  model 
sediment  with  variations  in  prototype. 

Grain  Size  (de),  mm 

Phi 

Prototype 

Model 

41.8 

29.5 

2.19 

45.3 

29.5 

2.48 

48.9 

29.5 

2.81 

41.8 

23.6 

1.84 

45.3 

23.6 

1.98 

48.9 

23.6 

2.25 

41.8 

35.4 

2.63 

45.3 

35.4 

2.97 

48.9 

35.4 

3.37 

41.8 

41.3 

3.07 

45.3 

41.3 

3.47 

48.9 

41.3 

3.94 

From  the  above  analysis,  a  very  fine  gravel  material  was  selected  as  the  best  option  for  the  1 :20 
model.  The  selected  material  has  a  D50  of  3.2  mm,  DIO  of  2.5  mm,  and  a  D90  of  4.6  mm.  Thus, 
the  lambda  (1),  or  scale  ratio,  for  the  models  BML  is  0.11.  This  1  is  much  greater  than  the  1 
selected  for  the  model,  but  this  is  not  of  concern  since  the  Shields  criterion  only  requires  that  the 
critical  shear  stress  is  scaled.  The  effective  roughness  height  for  prototype  is  based  on  the  D90 
(van  Rijn  1984).  The  average  D90  for  the  prototype  is  100  mm,  so  the  1  for  model  roughness  is 
0.046.  This  is  slightly  less  than  the  geometric  1  of  0.05  (1:20  scale)  but  is  arguably  sufficient  due 
to  the  physical  limitations  of  the  model. 

WASH  LOAD  SCALING:  The  WL  was  scaled  for  Rouse  similitude.  The  assumption  is  that 
there  is  enough  energy  in  the  system  (model)  to  keep  the  WL  suspended  until  it  encounters  an 
area  of  momentum  loss  (e.g.,  behind  the  barrier).  Satisfaction  of  Rouse  similitude  requires  that 
the  settling  velocity  is  scaled  such  that  Froude  similitude  is  satisfied.  This  will  also  ensure  that 
the  model  has  consistent  scaling  with  the  vertical  velocity  and  the  shear  velocity.  Froude 
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similitude  requires  that  the  velocity  scales  as  the  square  root  of  the  length  scale  (Equation  [8]). 
Hence,  Rouse  similitude  requires  that  the  settling  velocity  of  the  model  WL  sediment  scales  as 
the  square  root  of  the  settling  velocity  of  the  prototype  WL  sediment. 


w  =—. — w  „  =0.2236 w  „ 

a/  20  p'  p‘ 


(8) 


The  effective  grain  size  of  the  WL  in  the  prototype  is  0.26  mm  (settling  velocity  is  32.7  mm/s, 
Julien  [2002]).  Scaling  the  fall  velocity  of  the  effective  size  results  in  a  model  WL  grain  size  of 
0.1  mm.  However,  a  grain  size  of  0.1  mm  is  a  fine  sediment  and  is  challenging  to  procure.  The 
closest  available  material  has  a  D50  of  0.17  mm,  D10  of  0.11  mm,  and  a  D90  of  0.28  mm.  The 
settling  velocity  of  this  material  is  17.7  mm/s  compared  to  the  required  7.3  mm/s  (Julien  2002). 
Alternative  materials  (specific  gravity  <  2.65)  were  also  evaluated  but  would  be  cost  prohibited 
due  to  the  volume  of  material  required  for  model  testing. 

CONCLUSION:  The  nature  of  the  sediment  in  the  prototype  creates  an  advantage  for  sediment 
scaling.  The  bimodal  behavior  of  the  prototype  sediment  allows  for  the  use  of  two  different 
sediment  sizes  in  the  model,  thus  allowing  for  scaling  one  representative  size  for  Shields 
(initiation  of  motion)  and  one  for  Rouse  (settling).  Congruently,  the  large  prototype  material 
allows  for  sediment  scaling  without  model  or  material  distortions.  This  allows  for  the  realistic 
three-dimensional  prototype  behavior  in  the  model,  yielding  a  properly  scaled  analysis  of  scour 
and  sedimentation  around  the  modeled  structure. 


Validation  of  the  model  sediment  behavior  was  limited  to  qualitative  observation  from  site 
experts  (both  Seattle  District  hydraulic  engineers  and  the  fish  barrier  operators).  At  a  flow  of 
approximately  4,000  cubic  feet/second,  the  bed  mobilizes  and  readily  moves.  A  similar 
mobilization  occurred  in  the  model  for  the  BML  model  sediment.  While  not  definitive,  the 
similar  BML  mobilization  point,  between  prototype  and  model,  provides  some  confidence  in  the 
scaled  sediments  performance.  Also,  there  has  been  no  validation  for  the  WL  model  sediment. 

Application  of  the  method  should  be  limited  to  physical  modeling  of  granular  non-cohesive 
sediments  where  a  bimodal  prototype  gradation  exists.  Alternatively,  if  there  is  no  bimodal 
gradation  but  a  dominance  or  interest  in  studying  either  initiation  of  motion  or  settling,  then  the 
method  for  the  desired  behavior  can  be  used  independently. 

NOMENCLATURE 


d  =  grain  size 

de  =  effective  grain  size 

ds  =  grain  size 

d*  =  dimensionless  grain  size 

g  =  acceleration  of  gravity 

y  =  specific  weight  of  water 

ys  =  specific  weight  of  sediment 
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A  =  geometric  scale  ratio 

IRe*  =  particle  Reynolds  number 

SG  =  specific  gravity 

tc  =  critical  shear  stress 

9C  or  t*c  =  Shields  number 

[/*  =  friction  velocity 

v  =  kinematic  viscosity 

P  =  Rouse  number 

<p  =  Angle  of  Repose 

y  =  percentage  passing  for  ds 

ws  =  settling  velocity 

ADDITIONAL  INFORMATION:  This  technical  note  was  prepared  by  Jeremy  A.  Sharp, 

research  hydraulic  engineer  at  the  U.S.  Army  Engineer  Research  and  Development  Center, 

Coastal  and  Hydraulics  Laboratory.  Questions  about  this  technical  note  can  be  addressed  to  Mr. 

Sharp  at  601-634-4212  or  /ereffly.A£7tarp(fltoace.amv.ffl//. 

This  CHETN  should  be  cited  as  follows: 

Sharp,  J.  A.,  G.  L.  Brown,  and  G.  L.  Bell.  2017.  Sediment  scaling  for  Mud 
Mountain  Fish  barrier  structure.  ERDC/CHL  CHETN-VII-16.  Vicksburg,  MS: 

U.S.  Army  Engineer  Research  and  Development  Center.  http://dx. doi. org/ 
HL21079/1J_68 1/22 649. 
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